Diffraction of He atoms from Xe monolayer adsorbed on the graphite (0001) revisited: 

The importance of multiple scattering processes 
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We comment and discuss the findings and conclusions of a recent theoretical study of the diffrac- 
tion of He atoms from a monolayer of Xe atoms adsorbed on the graphite (0001) surface [Khokonov 
et al., Surf. Sci. 496(2002)L13]. By revisiting the problem we demonstrate that all main conclusions 
of Khokonov et al. that pertain to the studied system are at variance with the available experi- 
mental and theoretical evidence and the results of multiple scattering calculations presented in this 
comment. 
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In a recent Lettcr'^ and subsequent Erratum^ Khokonov, Kokov and Karamurzov^ (hereafter to be referred to as 
KKK) treated the problem of He atom scattering (HAS) from an ordered monolayer of Xe atoms adsorbed on the 
(0001) surface of graphite (Gr). The authors discussed a specific model of phonons in the Xe monolayer and used it 
in a scattering calculation based on the hard-wall model and the eikonal approximation with the aim to reproduce 
and interpret the experimental He atom diffraction intensities reported by Bracco et al.^. However, as we shall show, 
neither their description of phonons in the Xe overlayer, nor the results of calculations for the scattering intensities 
are consistent with the previous experimental and theoretical studies of the system Xe/Gr(0001) and the results of 
revisited calculations described below. 

The measurements of diffraction intensities from He beams incident normal to a (\/3 x \/3)R30^ monolayer lattice of 
Xe atoms on Gr(OOOl) carried out by Bracco et al."^ indicated large corrugation amplitudes as probed by HAS. These 
authors have attempted to interpret the measured diffraction spectra by using the hard corrugated wall (HCW) 
to model surface corrugation amplitudes obtained from pairwise summation of atomic He-Xc potentials^, and the 
eikonal approximation^ to calculate the diffraction intensities pertinent to the thus constructed HCW. However, their 
calculations showed that a more realistic scattering model was needed to obtain a better agreement between the 
experimental data and theoretical results. 

The studies of phonons in a Xe monolayer on Gr(OOOl) surface by HAS have been presented in Ref. 6. These 
measurements have demonstrated that the scattered He atoms couple most strongly to a low energy dispersionless or 
Einstein-like mode with frequency of ~ 3.2 meV. Both the lattice dynamical calculations'" and the molecular dynamics 
simulations^ utilising Xe-Xe potentials known from the gas-phase have shown that this nearly dispersionless mode 
corresponds to the vibrations of Xe atoms that are polarised perpendicular (vertical) to the surface plane, i.e. to a so 
called FT^ or S-phonon mode localised in the adlayer (for illustration of dispersion of all three Xe adlayer modes see 
Fig. 4 of Ref. 8). Low energy modes of the same S-character have been observed and measured also in Xe monolayers 
on Cu*^, Pt^°, Ag^^ and NaCl^^ surfaces. It has also been shown^'^ that a mode of this character hybridises with 
the substrate Rayleigh wave mode only for the two-dimensional phonon wavevectors restricted to a small region 
around the avoided crossing of the two dispersion curves in the first surface Brillouin zone (SBZ) where the S-mode 
polarisation vector is no more strictly localised in and perpendicular to the adlayer. 

In their treatment of Xe overlayer phonons, KKK have assumed a rigid Gr substrate and in this approximation 
the S-mode remains localised in the adlayer for all phonon wavevectors in the; SBZ. Their calculations predict that in 
Xe/Gr(0001) the vertically polarised S-mode exhibits dispersion from the value of 0.8 meV at the centre of the first 
SBZ up to 4.5 meV at the zone edge (see Fig 2. of Ref. 1). A molecular dynamics simulation of the same system, that 
was also based on the assumption of a rigid substrate, was presented in Ref. 8 and yielded an almost dispersionless S- 
mode with vertical polarization, in accord with the experimental evidence^. Lattice dynamics calculations of similar 
Xe adlayer systems, i.e. Xe/Cu(100) and Xe/Cu(lll), have been presented in Ref. 9. There, it was found first 
experimentally, and then confirmed theoretically that also in these systems the adlayer S-mode is to an excellent 
approximation completely dispersionless. Hence, all the results pertaining to vibrations of Xe adlayers on Gr(OOOl) 
and several other substrates are at strong variance with the first finding of KKK who obtain that the S-mode in 
Xe/Gr(0001) exhibits strong dispersion. This type of acoustic-like dispersion is typical of the adlayer modes with 
longitudinal (L) and shear horizontal (SH) polarisation that couple much weaker to the scattered He atoms^^ and 
hence their presence cannot explain why in the experiments of inelastic HAS from Xe/Gr(0001) surface^ the coupling is 
strongest to an almost dispersionless mode with a frequency of 3.2 meV. Hence, as a consequence of the inappropriate 
description of Xe adlayer in Ref. 1, the density of states they obtain for the S-mode displays a Debye-like behaviour 
(Fig 3. a of Ref. 1), whereas it should follow the density of states characteristic of a fiat, dispersionless mode. Thus, 
the mean square amplitudes of Xe atoms KKK could derive from such an inadequate density of states (Fig 3.b of Ref. 
1) would be generally wrong, although their magnitudes in a particular case may be close to the correct ones due to 
the fact that the very different models of vibrations can produce similar mean square displacements. 

Although the above described erroneous premises and lattice dynamics description of Xe/Gr(GG01) employed by 
KKK suffice to disqualify all the main conclusions of their analysis^, we also show that their scattering calculation is 
inadequate and hence cannot provide a realistic description of the collision system He^Xe/Gr(0001). 

KKK have used^ an eikonal-like semiclassical expression^ to calculate the scattering amplitudes in He^Xe/Gr(0001) 
collisions and in this approach the Xc/Gr(0001) surface was modelled by a corrugated hard wall without introducing 
Beeby's correction to partly compensate the neglect of the potential well in front of the surface. However, the 
applicability of this approach suffers from two fundamental restrictions^^ which were already noted and discussed in 
Ref. 3. First, the eikonal approximation is known to neglect the multiple scattering processes and therefore cannot 
be used in the cases of large surface corrugations and low projectile incident energies in which multiple scattering 
becomes important. Second, the calculations based on the HCW model can serve only as a rough estimate of the 
corrugation amplitudes and cannot be used as a tool for quantitative determinations of the surface structure and 
the details of interaction potential. These restrictions become especially important in the case of strongly corrugated 
surface profiles as is the case with the present system (see Fig. 1). In particular, if the repulsive part of He-target 
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potential is treated as an infinitely steep wall it is known that it cannot realistically represent the features of He-target 
interactions^^. Also, since in the HCW model the attractive component of the interaction is completely neglected, 
this approximation cannot be expected to hold at low incident He atom energies. Even for the room temperature 
He beam energies (the case considered by KKK), for which £"^=63. 8 meV, this approximation cannot be expected 
to provide reliable values for the intensities of diffraction peaks corresponding to He atoms scattered into the states 
with low energy in the 2;-direction (perpendicular to the surface). The well depth D of the He-Xe/graphite potential 
is known with high precision from the studies reported in Rcfs. 17 and 18 and amounts to D = 7.4 meV, and hence 
the He atoms scattered into the (50) and (60) channels (c.f. Figs. 1 and 4 of Ref. 1 and Fig. 2 of Ref. 3) have 
components of energy in the z-direction of 26 meV and 9.4 meV, respectively, that are comparable to D. To remedy 
the shortcomings of the HCW model in applications to the He^Xe/Gr(0001) system, KKK have assumed^'^ the 
magnitudes of corrugation parameters of the HCW shape function ((R) for modelling the Xe/Gr(0001) surface that 
are somewhat different from the best fit parameters suggested by Bracco et al.^ and, moreover, are inconsistent with 
the corrugation of the isopotential surface at 64 meV that can be constructed from the best available He-Xe gas phase 
binary potentials (c.f Fig. 1). Also, the HCW shape function employed by KKK^'^ in their eikonal approximation 
calculation of the scattering intensities, despite having the same functional form as the one given by Bracco et al., 
may appear as a possible source of computational errors because it is manifestly expressed in terms of the rectangular 
rather than oblique Cartesian coordinates introduced in Ref. 3 (c.f. Fig. 1 of Ref. 1 and inset in Fig. 2 of Ref. 3). 

To quantitatively substantiate the inadequacy of the scattering model on which the calculations in Refs. 1 and 2 
had been based, we have performed a coupled channel calculation^^ (CC) of the diffraction intensities for helium atom 
scattering from Xe/Gr(0001) surface^^ using the realistic He-target potentials known from the literature. We treat the 
target as static and construct the total He-Xe/Gr(0001) potential as a sum of the pairwise He-Xe interactions known 
from the gas phase^° and the long range interaction of He with the graphite substrate^^'^^. This potential is very 
similar to the one used by Hutson and Schwartz-'^^, and the very small differences are due to the fact that the He-Xe 
gas-phase potential we use is the one suggested by Cvetko et al.^", whereas Hutson and Schwartz^^ used a slightly 
different gas-phase He-Xc potential. The obtained total potential gives rise to the equipotential surface at 64 meV 
which exhibits the peak-to-peak corrugation amplitudes xi = 0.74 A and X2 = 0.96 A in the two high symmetry 
directions along the surface (see inset in Fig. 1). A summary of peak-to-peak corrugation amplitudes characteristic 
of the HCW models of Refs. 1 and 2, of Ref. 3, and of the ones obtained from the present calculation is given in 
Table 1. Note here that in order to obtain the amplitudes given in the first row of this Table we have corrected two 
errors appearing in the HCW shape function Co(R-) quoted in Eq. (7) of Ref. 1, viz. the misprinted value of the 
parameter Cia taken from Ref. 3, which should read Cia = 0.098 A(c.f. Ref. 2), and the role of x and y coordinates in 
the shape function of Ref. 1 which were here taken to have the same meaning of oblique coordinates as in Ref. 3 (c.f. 
the discussion at the end of the preceding paragraph). 

The CC calculations with a static periodic potential yield intensities of the (5-function-like diffraction peaks and 
in order to facilitate comparison with experiments we have broadened the calculated (5-functions by Gaussians of 
width aef to account for the finite energy spread of the incident beam, again given by a Gaussian of width aE- It is 
straightforward to show that the above Gaussian width parameters are related to each other, and for normal projectile 
incidence this relation reads: 

"^^^ ~ 2cos^/ ki Ei' ^ ' 

where G/ is the inverse lattice vector associated with a transition into a particular diffraction channel, ki is the 
projectile initial wave vector {ki = ^2mEi/h, m is the projectile mass) and 9f is the final scattering angle with 
respect to the ^;-axis (in our calculations, we have taken as/Ei = 7 %). 

In Fig. 2 we present the results of our CC calculations (full thick line in the lower panel) in comparison with the 
experimental data of Bracco et al."^ (upper panel). As can be seen from Fig. 2, the agreement with the experimental 
data is very satisfactory as regards the relative intensity ratios. The calculated relative intensities are also in a very 
good agreement with the results of earlier calculations by Hutson and Schwartz^^. To illustrate the importance of the 
attractive component of the interaction that was completely ignored in the calculations of KKK, we also display the 
results of independent CC calculations in which either of the two important features associated with the existence of 
the potential well and modifying the multiple scattering effects is neglected, i.e. (i) the exclusion of closed channels 
(evanescent waves) from the CC basis, or (ii) the exclusion of attractive component of the total potential from the 
calculation. As can be seen from the comparison of the full and dash-dotted line in lower panel of Fig. 2, the first effect 
caused by the exclusion of closed channels is large for the (60) diffraction peak whose intensity turns out to be more 
than a factor of 2 weaker in this approximation. However, the second effect of exclusion of the attractive potential is 
much more important since it influences the open diffraction channels in a more profound way and not only through 
the closed scattering channels. This effect is demonstrated by carrying out the CC calculations in which all the Fourier 
components of the attractive He-Xe/Gr(0001) interaction potential are neglected. The resulting diffraction spectrum 
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shown by the dashed hne in lower panel of Fig. 2 is very different from the one obtained with the full potential and 
does not reproduce the experimental spectrum in any of the important aspects. This is in full accord with Ref. 17 
where it has been shown that even small variations of the scattering potential can give rise to very large variations 
of the diffraction intensities. Hence, the recalculated CC-intensities that make use of the otherwise realistic He-Xe 
potentials with the attractive component removed produce the results which strongly differ from the ones obtained 
with the complete potential. This clearly demonstrates the inadequacy of the various approximate schemes based 
on the neglect of the effects of attractive components of the full potential in applications to He atom diffraction in 
the scattering regime in which the measurements reported in Ref. 3 were carried out. In other words, the relatively 
moderate differences among peak-to-pcak corugation amplitudes displayed in Table 1 may be deceptive as regards the 
behaviour of calculated diffraction intensities because they may give rise to different effects and results in different 
scattering models. 

The overall intensities of diffraction peaks obtained by the CC calculation with the full potential are about a factor 
of 10 larger than those obtained experimentally. This suggests that about 90% of He atoms in the experiment of Bracco 
et al.^ were scattered inelastically. Indeed, the experimental intensities when summed over all the measured diffraction 
peaks (c.f. Table 1 of Ref. 3) yield only 11.2 % of the incident intensity, which means that nearly 88.8 % of the scattered 
He atoms end up in inelastic channels, in excellent agreement with the factor of 10 that can be deduced from comparison 
of the results of our calculations and experiments. Note that this also implies a multiphonon scattering regime since the 
mean number of phonons n excited in a scattering event can be estimated from riexp = — lii(///o) = — In(O.ll) = 2.21, 
where I and Iq are the experimental values of the total intensity scattered into the elastic channels and the incident 
beam intensity, respectively, and I/Iq is the Debye- Waller factor (c.f. Ref. 13). It should be pointed out that this 
finding is again at variance with the conclusion of Ref. 1, that "the multiphonon processes can be neglected" . Namely, 
KKK have arrived at such a conclusion by uncritically employing Weare's criterion^^ for n in the regime in which this 
criterion is inapplicable (c.f. Figs. 1 and 2 in Ref. 22) and, moreover, made a numerical error in their estimate. To 
clarify this issue we have also calculated the theoretical value nth for the present system using the EBA multiphonon 
scattering formalism described in detail in Ref. 13. This yields rith = 2.6, again in a very good agreement with the 
above discussed value of riexp that can be extracted from experiments. 

In summary, by making use of the available He-Xe atomic pair potentials (Ref. 20) and the results of experimental 
investigations of the Xe/Gr(0001) surface by HAS (Refs. 3 and 6), the coupled channel method^^ for calculation 
of the diffraction intensities in thermal energy atom scattering from corrugated surfaces, and of the multiphonon 
scattering formalism^'^^ in combination with the lattice and molecular dynamics analyses^'^ to calculate the Debye- 
Waller attenuation of diffraction intensities, we have demonstrated that the results and conclusions reached by the 
authors of Ref. 1 and pertaining to: (i) the dispersion of vertically polarised modes in the Xe monolayer, (ii) the 
properties and modelling of the He-Xe/Gr(0001) interaction with the relatively large potential well and surface 
corrugation profile, and (iii) the role of multiple scattering processes in thermal energy HAS from Xc/Gr(0001) 
surface, are cither inadequate or in a complete disagreement both with the existing experimental evidence and with 
the theoretical assessments of the He^Xe/Gr(0001) collision system. The arguments which we supply to refute the 
conclusions reached in Ref. 1 show that: (i) the vertically polarised vibrational mode in Xe monolayer on Gr(OOOl) 
surface is nearly dispersionless, i.e. it is a flat S-mode which due to its low excitation frequency and perpendicular 
to the surface polarisation strongly couples to the scattered He atoms, (ii) the potential well depth and the surface 
corrugation amplitudes characteristic of the He-Xe/ Gr (0001) interaction under the studied scattering conditions are 
large, which necessitates a self-consistent treatment of the diffraction intensities within a multiple scattering formalism 
(and therefore any agreement between the results of the diffraction intensity calculations of Ref. 1 and experimental 
data can be only accidental), and (iii) the multiphonon processes arising from the strong He atom-S-mode coupling 
play an important role in the Debye- Waller factor-induced attenuation of experimental diffraction peak intensities 
reported in Ref. 3 with which the authors of Refs. 1 and 2 have compared their calculated intensities. 



1 A. Kh. Khokonov, Z.A. Kokov and B.S. Karamurzov, Surf. Sci. 496(2002)L13. 

^ A. Kh. Khokonov, Z.A. Kokov and B.S. Karamurzov, (Erratum to Ref. 1), Surf. Sci. in press. 
^ G. Bracco, P. Cantini, A. Glachant and R. Tatarek, Surf. Sci. 125(1983)L81. 
* G. Scoles, Ann. Rev. Phys. Chem. 31(1980)81. 

^ U. Garibaldi, A.C. Levi, R. Spadacini and G.E. Tommei, Surf. Sci. 48(1975)669. 

® J.P. Toennies and R. VoUmer, Phys. Rev. B 40(1989)3495. 

^ E. De Rouffignac, G.P. Alldredge, F.W. De Wette, Phys. Rev. B 24(1981)6050. 



4 



* M. Marchese, G. Jacucci and M.L. Klein, Surf. Sci. 145(1984)364. 

^ A. Siber, B. Gumhalter, J. Braun, A.P. Graham, M.F. Bcrtino, J. P. Toennies, D. Fuhrmann, Ch, W611, Phys. Rev. B 
59(1999)5898. 

° L.W. Bruch, A.P. Graham and J.P. Toennies, J. Chem. Phys. 112(2000)3314. 
1 K.D. Gibson and S.J. Sibener, J. Chem. Phys. 88(1988)7862 

^ R. Gcrlach, A.P. Graham, J.P. Toennies, H. Weiss, J. Chem. Phys. 109(1998)5319. 
^ B. Gumhalter, Physics Reports 351(2001)1. 

D. Farias and K.-H. Rieder, Rep. Prog. Phys. 61(1998)1575, and references therein 
^ G. Vidah, G. Ihm, H.-Y. Kim and M.W. Cole, Surface Science Reports 12(1991)133. 
^ A. Siber, Ph.D Thesis, University of Zagreb, 2002, unpublished. Downloadable at 

http : //bobi . if s .hr/~asiber/PhDThesis/ 

J.M. Hutson and C. Schwartz, J. Chem. Phys. 79(1983)5179. 
^ R.A. Aziz, U. Buck, H.Jonsson, J.-C. Ruiz-Suarez, B. Schmidt, G. Scoles, M.J. Slaman and J. Xu, J. Chem. Phys. 
91(1989)6477. 

^ G. Wolken Jr., J. Chem. Phys. 58(1973)3047; ibid. 59(1973)1159. 

'■° D. Cvctko, A. Lausi, A. Morganto, F. Tommasini, P. Cortona and M.G. Dondi, J. Chem. Phys. 100(1994)2052. 

'■^ J.H. Weare, J. Chem. Phys. 61(1974)2900. 

'■^ A. Siber and B. Gumhalter, Surf. Sci. 385(1997)270. 



FIG. 1. Contour plot of the Ho-Xc/Gr(0001) interaction potential y(R, z) (in rncV) along two high surface symmetry lines 
(a) and (b) shown in the inset. The energies of the isopotential lines are given in units of meV. The thick contour denotes the 
isopotential surface for which V{R, z) = Ei = 64 meV. 



FIG. 2. Top panel: Experimental diffraction intensities for He atom scattering from Xc/Gr(0001) surface shown as a function 
of the final scattering angle (after Ref. 3). Scattering parameters given in the inset. Bottom panel: Full line: Results of converged 
CC calculations using the total He-Xe/Gr(0001) potential as described in the text. Dash-dotted line: Results of CC calculations 
excluding closed diffraction channels. Dotted line: Results of CC calculations using the He-Xe/Gr(0001) potential without the 
attractive components. 



Reference 


XI [A] 


X2 [A] 


KKK'-^ 


0.72 


0.88 


Bracco et al.'^ 


0.64 


0.97 


This work 


0.74 


0.96 



TABLE 1. First and second row: pcak-to-poak HCW corrugation amplitudes used to compute He beam diffraction intensities 
in Refs. 1 and 2, and in Ref. 3, respectively. Third row: the corresponding amplitudes derived in the present work from the 
classical turning points in the He-Xe/Gr(0001) interaction potential V^(R, z) for He atom incident energy of 64 meV (see Fig. 
1)- 
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